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Abstract 

Three pion dominated observables of the parity nonconserving interactions between the cold neu- 
trons and parahydrogen are calculated. The transversely polarized neutron spin rotation, unpolar- 
ized neutron longitudinal polarization, and photon-asymmetry of the radiative polarized neutron 
capture are considered. For the numerical evaluation of the observables, the strong interactions 
are taken into account by the Reid93 potential and the parity nonconserving interactions by the 
DDH model along with the two-pion exchange. 
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I. INTRODUCTION 



The knowledge of the strangeness conserving hadronic weak interaction rehes completely 
on the parity nonconserving (PNC) observables. The PNC two-nucleon (NN) interactions 
provide a possible access to this least understood sector of the Standard Model. Experi- 
mentally such subtle particle spin control based PNC measurements are feasible but highly 
demanding. On the theoretical side, the challenge lies largely in the poorly known coupling 
values which parametrize the strength of the minuscule-sized PNC signal. 

The PNC A^A^ interaction is compelled to change either the spin or isospin of the system 
due to the Pauli exclusion principle and consequently the potential of the interaction is com- 
posed of various spin-isospin operators weighted by coupling constants. There are at least 
three considerable candidates for the PNC A^A^ potentials. The new model-independent 
effective field theory (EFT) approach offers two alternative choices for these potentials, 
namely the pionless and pionful The pionless one comprises only the short-range 

contact interaction whereas the pionful also the long- and medium-range interactions me- 
diated respectively by the one- and two-pion exchanges. The third potential is the most 
conventional DDH meson-exchange model 3| which takes into account the long- and short- 
range effects in terms of the single 7r^,p, and u exchanges but not the two-pion exchange 
contributions which are supposedly important in the medium-range. In spite of the different 
approaches, the operators appearing in the potentials are essentially the same, except that 
the ranges of the force in them varies. Another similarity with the potentials is that they are 
all parametrized in terms of about half a dozen ill-known couplings. In EFT these so-called 
low energy constants are expected to be extracted from the experimental data of a series 
of prospective high-precision measurements. As a downside, the EFT potentials cannot be 
used in the evaluation of observables yet. Therfore, for over thirty years up until today, the 
theoretical predictions have largely rest on the DDH model and their recommended "best" 
values for the weak couplings. The PNC one-pion exchange potentials of the pionful EFT 
and DDH model coincide and are proportional to the weak NNtt coupling The strength 
of the PNC two-pion exchange is also dependent on the same coupling. Besides the DDH, 
there are various calculations j^-H] for the the h^'^ (ranging between and 3.4 x 10""^) 



indicating a smaller value than what is the DDH "best" recommendation. The ongoing 



NPDGamma experiment 



12| is hoped to reduce the vagueness related to the h^^ by mea- 



suring with high accuracy the 7-asymmetry of the np — )■ 'yd at threshold known to be nearly 
a 100% pion exchange dominated. 

Up till recently, the effect of the two-pion exchange has been considered small and ne- 
glected from the analyses of the PNC observables. However, recent calculations of the PNC 
longitudinal analyzing power Al in the pp elastic scattering show the importance of the 



two-pion exchange [l3|, The PNC pp reaction offers an auspicious opportunity to study 



the two-pion exchange contribution for a couple of reasons. Firstly, the single pion-exchange 
does not appear according to Barton's theorem 15|, from which it follows that the two-pion 
exchange represents the longest-ranged and probably dominated contribution. Secondly, 
there exist three high-precision measurements of the Al at different energies (Bonn at 13.6 
MeV, PSI at 45 MeV, and TRIUMF at 221.3 MeV) which can be compared to the theoret- 
ical predictions. Consequently, when it comes to the one-pion exchange dominant PNC np 
interactions, the effect of the two-pion exchange should also be counted in. 

When slow neutrons collide with hydrogen molecules, they either elastically scatter off 
or get absorbed in the protons resulting in deuterons and photons. This work presents 
calculations of the three different pion sensitive PNC observables arising from the cold 
neutron interaction with parahydrogen. Two are due to the PNC elastic scattering enabling 
the spin rotation ^0 and polarization of the neutrons and the third one is the 7- 
asymmetry in the radiative PNC capture of polarized neutrons. The spin rotation and 
polarization of the neutron in the PNC np scattering were first discussed in Refs. 
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the wavefunction based calculations have later been performed in Refs. 18142 1|. There has 
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also been experimental interest in measuring the PNC neutron spin rotation ^0 in a liquid 
parahydrogen target at the Neutron Spallation Source (SNS), see Ref. The radiative 



21 



23-I37 



PNC reaction np — ^7 is also discussed in multiple papers, of which Refs 
present numerical predictions. 

n 

In the present work the Reid93 potential [38'] is chosen to take care of the strong np 
interactions. The long- and medium-range parts of the weak interaction are respectively 
due to the one- and two- pion exchanges. The contributions of the two-pion exchanges are 
taken into account in the observables by using the two separate PNC NN two-pion exchange 
potentials taken from Refs. 



39 



40|. The former of these potentials counts in the A(1232) 
isobar effects while the latter one does not. The neutron spin rotation -^0 and polarization 



are given in terms of the reduced matrix elements of the spin-space operators, which 



are the basic building blocks of the PNC potentials. The radial Yukawa functions of the 
operators are used in unregularized form and parameterized in three different ways by the 
meson masses vr, p, and u. The calculations are performed in the distorted- wave Born 
approximation (DWBA), in which the phenomenological strong interaction wavefunctions 
sandwich the PNC operator. In the calculation of the 7-asymmetry the pion-exchange 
current effects are included in the form of two-body dipole operators. The required bound 
and continuum radial wavefunctions together with their tiny parity admixed components 
are obtained from the exact solution of the coupled Schrodinger equation. 

The remainder of the paper is organized as follows. Section [TTl gives the basic formalism 
of the cold neutron interaction with parahydrogen, Sec. IIIII presents the results of the 
calculated observables, and Sec. |IV] summarizes the work. 

II. FORMALISM 

Hydrogen exists in nature in molecular form, composed of two protons bound by two 
electrons. The hydrogen molecule (-^2) comes in two species called parahydrogen (iff) 
and orthohydrogen {H2) with proton spins aligned antiparallel and parallel respectively. 
The de Broglie wavelength of a neutron at energies of a few meV is much greater than 
the internuclear separation Rq = 0.75 A between the protons in the hydrogen molecule. 
Consequently the neutrons, if not captured by protons, scatter coherently off the two protons 
in the molecules. Parahydrogen molecule is spinless, since its proton spins couple to zero 
and, therefore, it cannot depolarize a polarized neutron when they scatter elastically. The 
protonic wavefunction of the hydrogen molecule must be antisymmetric from which follows 
that in the ground states of the para- and orthohydrogen molecules, the rotational energies 
are respectively zero and = 14.7 meV, where / = /ii?Q is the moment of inertia and n 
the reduced mass of two protons. This energy determines the upper limit of the neutron 
center of mass (CM.) energy in order not to get depolarized by the conversion of the para- 
to orthohydrogen molecule. 

The low energy neutron-parahydrogen {nH"^) interaction is a three-body problem that 
can crudely be simplified to a two-body neutron-proton (np) interaction problem. The np 
continuum wavefunctions, in which the z-axis is taken along the direction of fc, are of the 
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form 

{r\kz; Qmr,)^^^ = ^t^^f^L^ V (\mJfn^\SMs) 
kr ^-^ 

L'mp 

x(L05M5|JM^)41)(fc,r)^/;,^(f)(-)^+i|T0), (1) 

where the superscripts (±) refer to the incoming (— ) and outgoing (+) wave boundary 
conditions, "S^jifgir) are the eigenfunctions of the coupled angular momentum, and the 
quantum numbers LSJT are abbreviated to Q. The quantum numbers ST J do not change 
under strong interaction. Because of the antisymmetricity requirement of the wavefunction, 
the isospin T may as well be considered uniquely defined by the LS, and so, on occasion, 
the Q is also designated for convenience with the spectroscopic notation Q = An 
adequate deuteron wavefunction, in the 7- asymmetry A^, is composed of three partial waves, 
which are the usual tensor coupled ^5*1 — ^Di and tiny parity admixed ^Pi. The np bound 
state wavefunction is given as 

(r|M,) = J2 '^^^iX(^)m^ (2) 
with the normalization dr l^^iP = 1 and energy eigenvalue of —2.2246 MeV. 



A. Neutron scattering 

Since the interest here is in the coherent niff scattering, one must carefully take into 
account the relative motion between the neutron and system of the chemically bound pro- 
tons. When a low energy neutron comes across the molecules in the medium, it interacts 
collectively with a number of them. As a result, the scattered waves originating from the 
molecules, interfere with the through passing neutron and change its momentum. Applying 
the Lippmann-Schwinger equation for multiple point-like scatterers each located at Vj, the 
wave of a slow neutron after travelling through the target in the z-direction then becomes 

2 g I ■p — ■ I 

^iq'z ^ ^iqz ^ f^^^ J2 J '-e'^'^ , (3) 

where f{q,6 = 0) denotes the forward ni^l" scattering amplitude and q is the relative mo- 
mentum of the neutron and molecule. The sum of the spherical waves from the scatterers in 
Eq. ([3]) may be written in the form of an integral over a smooth distribution of scattering 
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centers in a cylindrical shaped target of infinite radius. For the neutrons travelling along 
the axis of the target, it then follows that the right hand side of Eq. ([3]) becomes e**^'^^~'^\ 
with 

^(g,^) = -^/(g,o), (4) 

where M is the particle density of the medium. The neutron thus gains the (p = {q — q')z 
amount of phase when propagating through a medium of length z. Equation (jl]) is related 
to the index of refraction n = q' / q = 1 — (p / qz m neutron optics. 

By the initial choice of the transversely polarized spin in the positive x-direction (a^.) = 
+1, the neutron spin wavefunction = (|+) + \—))/\/2 contains equal amount of ± 

helicities in the direction of its propagation along the z-axis. The PNC interaction favors one 
helicity state slightly more than the other and thus depending on this state, the neutrons 
scatter a bit differently. The neutron wavefunction accumulates the amount phase 
(where m„ is the spin polarization of the incident neutron) labeled individually for each two 
states when passing through the target. It follows straightforwardly from the expectation 
value of the spin (<t) , that the neutron spin rotates in the xy-plane around the z-axis if the 
real value of the subtraction between the helicity states of Eq. (jlj) is non-zero 

2tiMz 
Q 

In the case of the unpolarized neutron beam, the neutrons gain some amount of longitudinal 
polarization due to the parity nonconservation when propagating through a medium. The 
incident beam intensity loss is given by dl±{q,z) = —f/a±{q)I±{q, z)dz from which, with 
the help of the optical theorem, the fractional polarization follows as the difference between 
the I+{q, z) and /-(g, z) divided by their sum 

2ttMz 
Q 



^iq,z) = -^Re(/+i(g,0) -/_.(g,0)). (5) 



P(g, z) ^ -^Im (/^ . (g, 0) - /_ i (g, 0)) . (6) 



The niff scattering amplitude may be written as / = —{jl/^)a in which the /i ~ M/2 
and fi ~ 2M/3 are respectively the neutron-proton and neutron-molecule reduced masses 



with the average nucleon mass M = 939 MeV and d = (as + 3aj)/2 (see e.g. Ref. j4l|) is the 
coherent scattering length expressed in terms of the np scattering lengths for the singlet ^5*0 
and triplet ^Si channels. That is a = — //2, where / is the np scattering amplitude. The 
niff and np scattering amplitudes and momenta are related as / = {jl/2ji)f and q = {fl/^)k, 
where k is the relative momentum of the neutron and proton. The relevant PNC part of 



the forward np scattering amplitude fm„{k,0 = 0) in distorted wave Born approximation 
(DWBA) is given by 

fmAKO) = -^(-)(A;i;mjf™C|A;f;m„)W, (7) 
Zn 

where the matrix elements are Hermitian. The S ^ P transitions are sufficient in the low 
energy PNC scattering and also equally important in both ways. Considering the lowest 
amplitudes, the derivative of the common factor in Eqs. © f = ReO{k) and §^ £P = 
lmO{k) (both per unit length) becomes 

0{k) = {k) - V2W:^l (k) - Wl-l ik)) , (8) 

where the subscript ^5*1 includes also its tensor coupled partner, the ^Di partial wave. The 
matrix elements of the PNC potential are further written in terms of the matrix elements 
of the operators, which appear in the PNC potentials of the DDH model and EFT, as 

a 

a 

^■^si(k) = ]^E(^oTx].:^4Txi-(^) +^OH,XT,,(^)), (11) 

a 

where the meson label a (= p, uj) is for the DDH model. In the case of EFT the a's and 
summation symbols are omitted. The reduced matrix elements of the spin-space operators 

^qIoiJ^) = ^"^(^^; Q'\\{^i&^2) ■ HV,yjr)]±||fci; Q)W, (12) 

with (0 = ±, x) are separated into the commutator [©]_ and anticommutator [©]+ elements. 
The constants where T denotes the total isospin in the initial and T' in the final state, 

include the matrix elements of the isospin operators and the other parameters associated 
with the potential, e.g. in the case of the DDH and pion, the constant is Co]\.]_ = g^h^'^ /V^. 

B. Neutron capture 



The thermal neutron capture cross-section on molecular hydrogen does not depend on the 
interference or binding effects of the protons in the molecule 42|] . It is therefore sufficient to 

7 



FIG. 1. The diagrams for the magnetic dipole moments considered in the calculation of a{np — )• jd). 
From left to right they are called impulse (imp), seagull (sea), pion- in- flight (fly), and delta (A) 
diagram. The wavy line is a photon, the solid line is a nucleon, the dashed line is a pion, and the 
bar is a A-isobar. 

simply calculate the neutron capture cross-section on free protons. The Ml ^Sq — ?■ ^Si — ^Di 
transition dominates the np — )■ '-yd reaction at threshold. By far the largest contribution (of 
about 90 %) of this reaction arises from the impulse approximation which couples the S- 
states. However, the one-pion exchange currents can also reach the D-state of the deuteron 
and play an important role in explaining the experimental value of the cross-section for 
thermal neutrons as was shown in Ref. 431]. 

The relevant photoproduction vertices, in terms of Lagrangian densities, are for the 
'jNNn interaction 



(13) 



and for the 'jtctc interaction 

= -e{d^77 X 77),A^. (14) 

Since the energy of the resulting photon at threshold of the reaction is only about 2 MeV, 
its wavelength is much larger than the deuteron size, and thus the electric E and magnetic 
B fields can be taken as constants. The scalar and vector potentials of the uniform (static) 
fields E and B are (plr) = —E ■ r and A{r) = ^B x r respectively. 

A diagrammatic illustration for the one- and two-body magnetic dipole moment operators 
is given in Fig. [H Besides Eqs. ( !T3|) and f fT4l) . the other necessary 'yNA, ttNA, and (PC and 
PNC) ttNN Lagrangians are given (in nonrelativistic form) in Appendix El In the impulse 
approximation, the relevant spin changing part of the operator is 



rn^"P = ^(ri.-f2,)(^i-^2). 



(15) 



The seagull contribution follows from Eqs. ( 1T3|) and (lAip leading to the exchange operator 



" 2 Vm, 



Tl X T2)^|f[f ■ (cTi X (T2)] -(TiX (T2^{1 + m^r) 



(16) 



Similarly, the pion-in-flight contribution (the Sachs exchange moment) is given by Eqs. (H 
and flAB 



rn^y(r) 



2 Vm 



) (ti X T2)z^f[r ■ (cTi X (T2)]{l + m^r) + cti x (T2{1 - m^r)|' 



(17) 



The final correction becomes by taking into account the static A^A intermediate state by 
using Eqs. (jSI]), dM]), and ( lA3ll 



m^fr) 



(ijiz - T2z) - i{ri X T2)^)|i(cri x f)(cr2 ■ f) 



9 (Ma - M) 
- i{a2 X f)(cri ■ f) - 2{ai - (T2) ■ rr 



(18) 



where the irrelevant isospin conserving terms are omitted as in the case of m'™P. The total 
magnetic moment operator is given by rh = Th™P + rh^^^ + rh^^ + rh^. 

The nonzero 7-asymmetry A-y arises from the interference between the Ml transition 
and the PNC interaction propelled El transitions. The El transitions connect the initial 
^Si — ^Di and final deuteron states through the parity admixed continuum and bound ^Pi 
states, while the PNC El transitions connect them directly. Figure [2] gives the diagrams for 
electric dipole moment contributions. The one-body operator is given by 



r2z)r 



and the two-body PNC exchange operator results from Eqs. (IT^ and (lA4p 



*2 firh-jT 



[Ti-T2-TizT2z 



AMm^V2 
The total El operator is the sum fie 



2(0-1 + 0-2) 



imp 



+ /i™^. Note that the changes parity, 



■r(o-i + o-2) 



(19) 



■ (20) 



while the /i™*" conserves it. Reference 4^ investigates the PNC El transitions and gives 
an additional 'jNNtt vertex leading to the spin-changing PNC El operator. However, this 
vertex has a vanishing contribution to the because the term proportional to the B is 
negligibly small and the Ml-El interference disappears as a consequence of the initial ^Sq 
states in both amplitudes. 

In terms of the reduced matrix elements of the electric and magnetic dipole transition 
operators, the 7-asymmetry A.y = {da^i — dcr_i)/{da_^i + da_i) reads now 

Ell^ J,- drV,. (r) (3Lf 1 1 Ae(r) | ^) " 



Ay{k) = V2Re 



^ El^ r d^'^L^i^) i'LiWfiUr) 1 1 i5o)wgj (r, k) 



(21) 
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FIG. 2. The El contributions to the 7-asymmetry A^. The black dot denotes the parity noncon- 
serving vertex. 

where fim = erh/2M and the multiplication factor cos 6', in which 6 is the angle between 
the neutron spin and photon direction, is left out. 



III. RESULTS 

The results for the neutron spin rotation -j^cj), polarization ^P, and 7-asymmetry 
are given in this section. The A^A^ wavefunctions are the exact solutions of the coupled 
Schrodinger equations where the strong distortions originate from the Reid93 potential 38 |. 
The OME parts of the observables can be and are, as customary, evaluated without regular- 
ization. The ^0 and are low energy scattering problems which are handled perturba- 
tively within the DWBA and therefore the inclusion of the form factors is not crucial. The 
DWBA approach allows the PNC potentials to be omitted from the Schrodinger equations 
because of their diminutive distortive effects to the strong wavefunctions. This, however, is 
not possible in the calculation of the Aj since the parity admixed wavefunctions are needed. 
In this case also the regularization of the PNC TPE potentials becomes necessity for the 
reason that they are too singular to be treated without it. As for the DDH model results, 
the DDH "best" values for the weak couplings are used along with the strong couplings 
Qt, = 13.45, Qp = 2.79, and = 8.37 and the anomalies Xp = 3.71 and Xui = —0.12. 



The employed TPE potentials are taken from Refs. [391] and 



40| and respectively ab- 



breviated as V^k"^^(*") and V^^^{r) of which the former takes into account the A^A^ and 
A^A intermediate states while the latter only the A^A^. The A^A^ intermediate parts of these 
potentials coincide, if the constant term in the dispersion relation of the V^^^{r) is ex- 
cluded. The constant term in momentum space corresponds to a term proportional to the 
delta function 6{r) in coordinate space. However, when the TPE potentials are regularized 
by form factors, their A^A^ parts inevitably differ from each other. The TPE potentials are 
modified by the monopole A^(qr^ + A^)~^ and dipole A^{q'^ + A^)^^ form factors and used 

10 



with three different cut-off mass values A = 0.8, 1.0, and 1.2 GeV. The same type of the 
monopole form factor emerged with the PNC TPE potential V^^^{r) was also used in Ref. 



34| to calculate the 7-asymmetry in -np — )■ 'yd at threshold. The results of the medium-range 



TPE are sensitive to the form factors, which have an increasingly suppressing effect on them 
when the cut-off is decreased and the rank of the form factor is raised. 



A. Neutron spin rotation -^(p and polarization -^P 

The spin-space operators appearing in Eq. f|T2l) are common building blocks in the PNC 
A^A^ interaction. In the approach in which the operators are placed between high-quality 
phenomenological wavef unctions, the low-energy < 1 keV scattering matrix elements can be 
expressed in constant form, as given in Tab. HI Therefore, by means of the matrix elements 
of Tab. m it is straightforward to customize the OME contributions within any model 
that uses the vr, u, and/or p exchanges as the ranges of the unmodified Yukawa functions 
y„(r) = e-™-747rr. 

The neutron spin rotation is split into one-meson exchange (OME) and two-pion exchange 
(TPE) components as ^0 = ■^<P'^^^ + ^0'^^^- Since a neutron scatters coherently from a 



TABLE I. The values of Eq. ([T2]) for the spin rotation components X|J"]_ (fe) = Re (A;"^ (fe)) 
and Xgj^j^(A;) = Iie(^—ik~^J'Q^^^^{k)^ in units of mb and for the spin polarization compo- 
nents ag^jJA;) = ^Im(fc-3j2'°_(A;)) and G^^^^Jk) = ^Im(-ifc-3 jg'"^ (A;)) in units of 

M 



10^1° fm/meV. The k and T£^^ (related as ^Laj^ = ^-) are respectively in units of fm ^ and meV. 



The functions I and G are constants within the neutron kinetic energy range of about T^sh~^'^ 
keV. 



5o[x 
So[- 

'3Sl[x 



-^Pia 



vr 



-30.1992 
-42.2133 
1.5588 
-4.4409 
3.4753 
5.6818 



P 



-0.7825 
-0.6494 
0.0868 
-0.1020 
0.1528 
0.0843 



-0.7444 
-0.6143 
0.0834 
-0.0975 
0.1459 
0.0798 





vr 


P 




r'-^Poa 


-8.6220 


-0.2234 


-0.2125 


^^Poa 


-12.0522 


-0.1854 


-0.1754 


Q^Pia 
^35i[x]_ 


-0.1020 


-0.0057 


-0.0055 




0.2905 


0.0067 


0.0064 




-0.2273 


-0.0100 


-0.0095 




-0.3717 


-0.0055 


-0.0052 
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pair of protons in the hydrogen molecule, it is then appropriate to use the molecule number 
density of liquid hydrogen instead of a two times larger atom number density. Therefore, 
the liquid parahydrogen particle density value of A/" = 0.021 molecules/A'^ is used in the 
numerical results. In terms of the DDH model, the rotation may be written as 

A^OME ^ /-Q.eir/iW - 0.138/ii°) - 0.012/i« - Q.l2Qhf + 0.004/i« + 0.130/if ) — . (22) 

This has the value of 3.31 x 10"^ when the DDH "best" values are plugged in. With 
this model and values, the one-pion exchange (OPE) has a dominance of about 85%. The 



result of Eq. 



is consistent with the one in Ref. 



20| , except half the size (because o: 



half the size particle density value), and also in line with the predictions of Refs. 



19 



the 



2l| 



which al 
in Ref. 



employ the Argonne fig potential. The result using the Paris potential reported 



18| is in the same order with the aforementioned results but of the opposite sign. 



Also in the result of Ref. jl8[ some of the signs between the partial contributions are in 



disagreement with the mutually consistent results of Ref. 



20| and Eq. (|22l). 



TABLE II. The TPE contributions to the neutron spin rotation -j^ff^^^ in units of /li^^^ and 
polarization Tf^4^P^^^ in units of /li^'' x 10^^^ fm^^m^^meV^^. The K and D stand for the 

-'Lab "-^ 

V^^^{r) and V^^^{r) while the F and FF stand respectively for the modification by K^{q^ + 
and A^(q^ + A^)~^, where the cut-off masses A are in units of GeV. 





A 


K(FF) 


D(FF) 


K(F) 


D(F) 


0.8 


-0.040 


0.026 


-0.105 


-0.020 


1.0 


-0.074 


0.004 


-0.132 


-0.040 


1.2 


-0.103 


-0.017 


-0.152 


-0.055 



k d pTPE 


A 


K(FF) 


D(FF) 


K(F) 


D(F) 


0.8 


2.648 


-1.712 


6.845 


1.295 


1.0 


4.844 


-2.442 


8.648 


2.619 


1.2 


6.763 


1.104 


9.925 


3.604 



Because the V^^^{r) and A^A^ part of the V^^^{r) are identical in unregularized form, 
they contribute the same amount to the rotation -j^4>Wn ~ — O.lOl/ii^'' ^ reducing the 
OPE effect by about 15%. By taking into account also the A effects of the V^^^(t'), the 
contribution doubles to j-^cfW^ = -0.209/il^^ '-^ cutting down the OPE effect by over 30%. 
However, since the TPE is a medium-range effect, it is sensitive to form factors. The form 
factor modified TPE contributions to the j^fp^^^ are given in Tab. Ull As is seen in the Table, 
the regularized TPE contributions are, like the unregularized ones, negative in almost all 
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cases and, thus, opposite to the OPE effect. Only the D(FF) case differs by having a different 
sign at small cut-off masses. 

When it comes to the neutron spin polarization ^P, the effect is much smaller than 
the one of ^0 and therefore not particularly interesting because it is experimentally much 
less achievable. For example, similarly to Eq. ( |22l) . the DDH OME contribution to the 
polarization with T^^b = 10 meV neutrons is only ^-PoDif — '''•43 x 10~^^ This result 
can be constructed by using the matrix elements ^q^qj^^I^) of Tab. [H which lead with the 
DDH "best" values to j^£pSdh = 8-17 x IQ-^^ fm^^m-^meV^^ The OPE contribution 
to this number in the same units is —18.54 and has, thus, about 70% dominance. As 
in the rotation, the unregularized TPE contributions to the j^-j^P'^^^ (in units of h^^ x 

Lab 

10^^^ fm~^m~^meV~^) are 6.63 and 13.70, in which the former is the same for both potentials 
when only the A^A^ intermediate states are considered and the latter is for the Vk"^^(^) when 
also the A^A intermediate states are included. As in the corresponding case of the rotation, 
the TPE effect with the A, diminish over 30% the OPE contribution. The effects of the 
form factors on the j^j^P'^^^ are presented in Tab. [Ill and show similar features as in the 

Lab 

case of the spin rotation. 



B. Photon asymmetry A^y 

An experimental value of the radiative thermal neutron capture cross-section of proton 
is (7{np — 'jd) = 334.2 ± 0.5 mb 45|. Theoretically it may be given by 

^cap = ^$^|y^ drV,,ir)eL,,\\ni{r)\\'So)Utlir,k) , (23) 



where a (= e ) is the fine structure constant, the CM. energy of the photon, and m(r) the 
magnetic moment operator, which is the sum of Eqs. fll5p - fll8l) . The used coupling values are 
= 13.45 iU = m^gj2M), = VW^U. and /ia = />./2^ {f^Ni^lm^ = enj2M), 
where g^^ is the irNN coupling, /* the quark model result for the ttNA coupling |46], /iA the 
transition magnetic moment, and fi^ = 4.71 the isovector magnetic moment of the nucleon. 
Evaluation of Eq. fl25]) for the thermal neutrons (T^^b ~ ^5 meV) gives (Xcap = 334.42 mb, 
which is in an excellent agreement with the experimental cross-section. In this result, the 
mere impulse approximation, i.e. rh(r) = rh™P, produces the largest part of the cross- 
section giving (T™^ = 303.81 mb. The missing ~ 10% enhancement results from the OPE 
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current corrections, as proposed in Ref. 43 1. 

By using the DDH model and their "best" couphng values, the OME contribution to 
the 7-asymmetry Eq. ( 121]) is ^^di^h = —5.387 x 10^^. Even though included, the effects 
of the PNC El current (~ — 0.2%o) and heavy meson p- and w-exchanges (less than 1%) 
are negligibly small. In terms of the weak pion coupling, the result may be written as 
-4°DDH ~ -0.117/ii^^ which is in harmony with the previous predictions (see the references 
mentioned in the Introduction) in most cases. The results, in which the TPE contributions 
are added on top of the ^^di^h; shown in Tab. [Till and are rather self-explaining. 
The result of the model D(F) with A = 1.0 GeV is consistent with Ref. 3J] and show 
~ 6% smaller 7-asymmetry than without the TPE. Just like in scattering cases, the TPE 
contribution of the D(FF) model with A = 0.8 and 1.0 GeV differs by the sign from the 
other models and, therefore, strengthens the total asymmetry. Otherwise the TPE effect 
is destructive as in the case of the ^0"""^^- The largest contribution arises from the K(F) 
model with A = 1.2 GeV, which diminishes the .A^j^!^ of about about 20%. 



TABLE III. Photon asymmetry A^y including the TPE in units of 10 



-8 





{np — )• 


7d) 




A 


K(FF) 


D(FF) 


K(F) 


D(F) 


0.8 


-5.08 


-5.58 


-4.60 


-5.24 


1.0 


-4.83 


-5.41 


-4.40 


-5.08 


1.2 


-4.61 


-5.26 


-4.26 


-4.97 



IV. SUMMARY 

Three PNC observables for cold neutron interaction with parahydrogen were calculated. 
All the observables, the neutron spin rotation ^0 and polarization in scattering, and 
the 7-asymmetry in capture, were found to be dominated by the pion exchange. The 
effect of the TPE was also taken into account in these observables and investigated in several 
settings. In all cases, the TPE effect was mainly opposite to the one of the OPE. 

The OME contribution to the ^0 was concluded to be a factor of two smaller than the 
most recent predictions and that the TPE decreased it up to ~ 30% further. The was 
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considered rather uninteresting due to its small size. In the A^, the OPE currents gave the 
expected increment for the Ml transitions but were insignificant for PNC El ones. The 
asymmetry was found to arise more or less completely from the pion exchange, i.e. the OPE 
weakened by the TPE up to ~ 20% or so. 

Unfortunately, so far there exist no direct and helpful experimental data for these specific 
observables. However, as already mentioned, there exist three precision experiment data 
points of the PNC longitudinal analyzing power Al{pp — > pp) of which two are low energy 
joints originating completely from the ^Sq — ^Po transition. A recent Al calculation of Ref. 



14j | uses exactly the same TPE models (potentials, couplings, and form factors) as used 
here. Based on the DDH model along with the TPE model K(FF) with A = 1.0 GeV gives 
the best fit and is also in a good agreement with the experimental data. This suggests that, 
in the calculated observables -^(p, ^P, and A-y, the destructive TPE effect compared to the 
OPE one is roughly 10% in each case. 
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Appendix A 

The nonrelativistic interaction Lagrangians (the couplings are explained in the text) are 

LnNN = — CT ■ Vt ■ TT, (Al) 

LyNA = ^S-Bt + R.c., (A2) 
L^NA = —S ■ VT ■ TT + H.C., (A3) 

LITn = ^(T X TT),. (A4) 

The S and T are respectively the NA spin and isospin transition operators [461] . 
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